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Abstract: In arXiv:1101.3326[hep-th], a (2+l)-dimensional holographic Josephson junc- 
tion was constructed, and it was shown that the DC Josephson current is proportional to 
the sine of the phase difference across the junction. In this paper, we extend this study 
to a holographic description for the (3+l)-dimensional holographic DC Josephson junc- 
tion. By solving numerically the coupled differential equations, we also obtain the familiar 
characteristics of Josephson junctions. 
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1 Introduction 

The study of the AdS/CFT correspondence and its applications have been the subject of 
great activity since the original conjecture was made by Maldacena [1]. The AdS/CFT 
correspondence is a powerful tool to study the strong coupled field theory. In the last few 
years, the application of the AdS/CFT correspondence to the Superconducting/Superfluid 
phase has been an interesting subject, which states that the Einstein-Maxwell-scalar model 
provides a holographically dual description of a superconductor [2-4]. For reviews on 
holographic superconductors, see [5-7]. 

Recently, in Ref. [8], the authors proposed a gravity model which provides a holo- 
graphically dual description of a 2+1 dimensional Josephson Junction. By choosing spatial 
dependence \i and non-vanishing constant J, which correspond to chemical potential and 
superfluid currents on the AdS boundary [9-12], respectively, the authors have solved 
numerically the coupled, nonlinear partial differential equations and reproduced the well- 
known result that the DC Josephson current is proportional to the sine of the phase dif- 
ference across the junction. Meanwhile, they also studied dependence of the maximum 
current on the temperature and size of the junction, which matches precisely with the 
results for condensed matter physics. 

In this paper, we would like to extend the work [8] to a holographic description for the 
(3+l)-dimensional holographic DC Josephson junction. We will study the 5D Einstein- 
Maxwell-scalar model and solve numerically the coupled, nonlinear partial differential 
equations of motion, and look for how to construct a holographic Josephson junctions 
by choosing the boundary conditions. 

The paper is organized as follows. In Sec. 2, we review 5D Einstein-Maxwell-scalar 
theories and set up a gravity dual of a (3+l)-dimensional SNS Josephson junction. In 
Sec. 3, we show numerical results of the EOMs and study the characteristics of the (3+1)- 
dimensional holographic Josephson junctions. The last section is devoted to conclusion. 
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2 The model 



In this section, we study the dynamics of a complex scalar field coupled to a U(l) vector 
potential in (4+l)-dimensional AdS spacetimes. The action of gravity can be written as 



S = J d 5 x^(R+^) , 



(2.1) 



where L is the curvature radius of asymptotic AdS Spacetimes. We are interested in the 
solution of the planar Schwarzschild black hole with the form: 



ds 2 



-f(r)dt 2 + 



dr 2 



+ r 2 (dx 2 + dy 2 + dz 2 ) 



M 



(2.2) 



where M is the mass of the black hole. The black hole has the event horizon rn at 

r H = M*L?. (2.3) 
The Hawking temperature of the black hole reads as 



T 



l_df_ 

4-7T dr 



r=r H 



M 1 / 4 
L 3 / 2 ' 



7T 



(2.4) 



In this background, we now consider a Maxwell field and a charged complex scalar field, 
with the action 



S 



d b x^g 



1 



F^F^ - \Vtp - iA^\ 2 - m 2 \4>\' 



where F„ v = du.A v — d v A^. The equations of motion include the scalar equation 

- (V - iA a ) (V a - iA a ) i> + mV = , 

and Maxwell's equations 

V a F ab = i [^*(V b - iA b )tp - V(V b + iA b )r\ ■ 

Taking a static ansatz 

<0 = \ip\e iip , A = A t dt + A T dr + A x dx , 



(2.5) 



(2.6) 



(2.7) 



(2- 



where <p, At, A r , and A x are real functions of r and x, one can introduce the gauge- 
invariant field M = A — dip, which is the same definition as that in [8]. 

With the black hole background (2.2) and the above ansatz (2.8), the equation of scalar 
field (2.6) can be written as the real part: 



1 



r 2 f 

and the imaginary part: 
1 



/ r 



1 / 'M? 



r2 Mi 



d r M r + -rrj.d x M x + — M r d r \4>\ + -£d. 



r 2 f 



f\f 



M x 



0, 



r 2 f 



(2.9) 



(2.10) 
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The imaginary part of the scalar equation can also come from the conservation of the source 
in Maxwell's equations. Meanwhile, Maxwell's equations (2.7) can be written as 

d 2 r M t + 4l%M t + -d r M t - m^ Mt = , (2.12) 
r J r J 

d 2 x M T - d r d x M x - 2r 2 \iP\ 2 M r = , (2.13) 

f 1 2h/;l 2 
d 2 M x - 8 r d x M r + ( J - + -) (d r M x - d x M r ) - -^M x = , (2.14) 
f r f 

where a prime denotes derivative with respect to r. Because Eqs. (2.9)-(2.14) are cou- 
pled nonlinear equations, one can not solve these equations analytically. However, it is 
straightforward to solve them numerically. 

In the following part of this paper, we choose the mass of the scalar field to be m 2 = —3, 
which is above the Breitenlohner-Freedman bound [14]. In order to solve these coupled 
equations, first, we need to impose regularity at the horizon and the boundary on the 
radial coordinate. At the horizon r = rjj, Mt(rn) = 0; At the Ads boundary, the scalar 
field takes the asymptotic form 

,^£M + ^ +0 (_L), ( , 15) 

the asymptotic behavior of the Maxwell fields are 

M t = Mx) -PM+ofy, (2.16) 

M r = (1) , (2.17) 

M x = v{x) + ^ + o[^ , (2.18) 

here, /i, p, v and J are interpreted as the chemical potential, charge density, superfluid 
velocity and current, respectively [13, 15]. J is a constant. For |^|, both tp^ and ip^ 
are normalizable, one can impose the condition either tjj^ or ip^ vanishes. For simplicity, 
we will consider the case tp^ = and interpret (O) = ip( 2 ^ as the condensate. Here, we 
introduce the same phase difference 7 = — j A x as the case in [8] , which is the gauge 
invariant and can be rewritten as 

/oo 
dx [u(x) - i/(±oo)] . (2.19) 
-00 

Second, on the spatial coordinate x, we impose the Dirichlet-like boundary condition 
on M r and Neumann-like boundary condition on Mt and M x at x = 0. At x = ±00, 
the function is ^-independent. Thus the boundary conditions of the coupled equations 
(2.9)-(2.14) are determined by J and /i. 

According to the study of paper [8] , we can also introduce the critical temperature T c 
of the junction , which is proportional to = /i(oo) = ^(—00): 



7T fl 



(2.20) 



c 
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where fj, c ~ 4.16. Furthermore, the effective critical temperature inside the gap is 

rp 1 MO) 
i = . 

7T [l c 



(2.21) 



In order to describe an SNS Josephen junction, we choose the below n(x) which is the same 
as that in [8]: 



1 - e 



2tanh(£) 



tanh 



x 



tanh 



a 



(2.22) 



where /ioo is the chemical potential at x = ±oo, and L is the width of our junction. The 
quantities a and e control the steepness and depth of our profile, respectively. 



3 Numerical results 

In this section, we will computer the solution of the coupled equations (2.9)-(2.14) numeri- 
cally. It is convenient to set the change of variables z = \ — rujr and x = tanh(^r). First, 
we show the solutions for and J\d x in Fig. 1, with /ioo 

= 6, L = 3, e = 0.6, and a = 0.5. 
Then, in Fig. 2, we find that the superfluid current is proportional to the sine of 
the phase difference across the junction, namely, the red dots come from the numerical 
calculations match precisely with the the black solid sine curve. Analyzing the graph, we 
can obtain the maximum current across the junction: 

J max /T 6 3 w 0.871. (3.1) 

In Fig. 3(a), the dependence of J max on the width of the gap is shown. The graph 
predicts an exponential decay with the growing width of the gap in J max : 

Jmax/T c 3 =A e~i. (3.2) 
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Figure 2. Superfluid current J max /T^ as the function of the phase difference 7. The black line is 
the sine curve. The parameters are set to /i,^ = 6, L = 3, e = 0.6, a = 0.5. 

In Fig. 3(b), the condensate (O) = tp^ at zero current is shown. The graph also predicts 
an exponential decay with the growing width of the gap in J max : 

(O) x=0 ,j= q /t3 = A ie ~^. (3.3) 

Fitting Eq. (3.2) and Eq. (3.3) with the two sets of data, we can obtain {£,Ao} ~ 
{0.85,30.32} and {£,^1} » {0.84,104.33} for Eq. (3.2) and Eq. (3.3), respectively. The 
disagreement of two values of £ is smaller than the 4D case in [8] . 

In Fig. 4, we obtain the relation of J max and T. Since e = 0.6, we show the region 
corresponds to T/T c < 0.6, which depicts the character of an SNS Josephson Junction. 

4 Conclusion 

In this paper, the holographic description for the holographic DC Josephson junction is 
extended to the case of the (3+l)-dimensional junction. In the (4+l)-dimensional gravity 
background, we investigate a set of five couple, partial differential equations. By choosing 
spatial dependence \i and non-vanishing constant J, we find the the superfluid current is 
proportional to the sine of the phase difference across the junction. Moreover, the graph 
which predicts an exponential decay with the growing width of gap in J max is obtained. 
At last, the relation of J max arid T is also obtained, from which one can verify that, near 
the critical temperature T c , J max can reach zero. 

From the result of this paper, one can see that the model of the holographic DC 
Josephson junction is valid for the (3+l)-dimensional case. It would be interesting to 
investigate the holographic DC Josephson junction in diverse dimensions and in various 
gravity theories. 
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Figure 3. The dependence of J max and (O) x= o on L. The parameters are set to = 6, e = 0.6, 
cr = 0.5. 
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Figure 4. The dependence of J max on T, The parameters are set to e = 0.6, a = 0.5 
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